Abstract 



In this paper we give definitions of basic concepts sucfi as symmetries, first integrals, Hamilto- 
nian and recursion operators suitable for ordinary differential equations on associative algebras, 
and in particular for matrix differential equations. We choose existence of hierarchies of first 
integrals and/or symmetries as a criterion for integrability and justify it by examples. Using 
our componentless approach we have solved a number of classification problems for integrable 
equations on free associative algebras. Also, in the simplest case, we have listed all possible 
Hamiltonian operators of low order. 
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1 Introduction 



In the classical theory (Lie, Liouville, etc.) of ordinary differential equations (ODEs) there are 
remarkable results which relate the property of integrability of ODEs in quadratures with the 
existence of continuous symmetries and first integrals. Symmetries and first integrals may serve 
as a solid mathematical foundation for an algebraic theory of integrable equations. In the case 
of integrable partial differential equations (PDEs) such a theory does already exist and proved to 
be extremely efficient ^ ^ ^ . The key property of integrable PDEs is the existence of infinite 
hierarchies of local infinitesimal symmetries generated by a recursion operator. Characteristic 
features of integrable Hamiltonian PDEs are multi-Hamiltonian structures and hierarchies of local 
conservation laws. These properties are well described in the fundamental monograph by P.Olver 
^ where references on original publications are well presented. 

A straightforward generalisation of these ideas to the case of ODEs is impossible, since the 
number of independent commuting symmetries and first integrals for a finite dimensional dynam- 
ical system is finite and bounded by its dimension. To overcome this obstacle we propose to study 
an intermediate object, namely, equations on free associative algebras. Here we are going to show 
that such equations are quite similar to PDEs. In particular they may have infinite hierarchies 
of first integrals or symmetries and there are plenty of reasons to choose these properties as an 
algebraic definition of integrability. 

If we are given a matrix differential equation with no restrictions on the matrix dimension 
then we can treat it as an equation on an abstract associative algebra. For example, equations 
for the classical Euler top in n-dimensional space can be written in a matrix form 

M = [M, Q\ . (1) 

Here M, are real skew-symmetric n x n matrices of angular momentum and angular velocity, 
brackets [ , ] denote usual matrix commutator and 

M = jn + nj, (2) 

where J is a constant diagonal matrix of moments of inertia. The problem of integration of the 
general n dimensional case was solved by Manakov [^. He made an important observation that 
the Euler equation (|l]), (|2|) is a stationary point of the A*'- wave equations, known to be integrable 

0- 

Equation (|^) is a Hamiltonian system^ with Hamiltonian 

Hi = itrace(Mf7) . (3) 



It is easy to check that 

ffa = -trace(M^J^), Hz = tra.ce{2J^ M'^ + f M .f M) (4) 

are also first integrals of (|^). First integral H2 generates a symmetry of the Euler equation 

Mr = [Af, gradAj-ff2] = J'^ -M'^J'^ . (5) 

In terms of M, J, Q, first integrals and symmetries have reasonably simple polynomial form. 
An attempt to find first integrals or symmetries in component notations, say by the method of 
indetermined coefficients, would fail for dimensions n > 4 because of enormous scale of computa- 
tions involved. In order to avoid component-wise computations and make results suitable for any 
dimensional problem we could regard objects AI, f2, J as generators of an associative algebra (i.e. 
assume that all our objects are polynomials of non-commutative variables M, fl, J with complex 
or real coefficients). This is not a free algebra because of the constraint (^). A more serious prob- 
lem is that the Euler equation (^) does not define a i-derivation of the algebra. Indeed, we cannot 
determine ^ as an element of the algebra using the Euler equation (|^ and the constraint 
In other words the Euler equation is not an evoJutionary equation. On the contrary, symmetry 
^ with the condition Jt = defines a derivation of a free algebra generated by elements M and 
J (we shall denote associative multiplication in this algebra by o). It is much easier to study (^ 



^Indeed, it can be rewritten in the form AI = {M,Hi} = [AI,gradj^jHi] with Poisson brackets defined as {F, 
trace (grad [M,gradjvfG]) . 



than the original Euler equation. Since the Euler equation and (^) belong to the same hierarchy, 
they have common symmetries and first integrals. 
In this paper we will use equation 

Mt = M'^ oC -Co (6) 

which coincides with (^) if C = J^, as a basic and instructive example. 

Being defined on a free algebra, equation has infinite hierarchies of polynomial symmetries 
and first integrals (see Section 1) and therefore is integrable according to our definition. Taking 
concrete finite and infinite dimensional associative algebras (such as algebras of differential or 
integral operators, matrix algebras or algebras with particular commutation relations) instead of 
the universal free associative algebra one can obtain a wide variety of equations which we expect 
to be integrable in a more conventional sense. For instance, in matrix realisations, symmetries 
and first integrals of equations on abstract associative algebras are inherited by the corresponding 
matrix ODEs and can be used for their integration in quadratures. 

Our approach to differential equations on associative algebras requires proper definitions for 
such concepts as symmetry, first integral, Hamiltonian and recursion operators, etc. In the next 
section we give a definition of symmetries and first integrals suitable for equations on free asso- 
ciative algebras. Also we show that equation (^) is not a unique representative of one component 
integrable equations and formulate a few classification results (we call them Results, rather than 
Proposition or Theorem, because they were obtained by straightforward computations). 

Frechet derivative, gradient, Hamiltonian and recursion operators for equations on free as- 
sociative algebras are defined in Section 3. In the same section we formulate a simple classi- 
fication result for Hamiltonian operators. In particular, we have shown that equation (|^ is a 
tri-Hamiltonian system. Ratios of the corresponding Hamiltonian operators yield two recursion 
operators. Existence of two independent recursion operators enables us to build up a complete 
two-index hierarchy of symmetries for equation (^). 

Many important multi-component integrable equations on associative algebras can be obtained 
as reductions of one-component equations. Here we give just one example. Let M and C in (^ 
be represented by matrices of the form 
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(7) 



where Uk and Jk are block matrices (of any dimension) or even generators of a bigger free algebra. 
It follows from equation (n) that satisfy the non-abelian Volterra equation 



— Uk = life O Uk + 1 O Jk + l 

at 



Jk-l O Uk~l O Uk, fc e Zjv 



(8) 



In Section 4 we consider other examples (modified Volterra equation, discrete Burgers type equa- 
tion, etc) and reductions to two-component systems of equations. Also we solve a simplest 
classification problem for two-component equations with one extra symmetry and study some of 
their properties. 

We conclude our paper by a discussion of possible developments and future research. 



1.1 First integrals and symmetries for equations on free associa- 
tive algebras 

Equation (^ is naturally defined on a free algebra M generated by elements M and C over the 
field C with an associative (but noncommutative) multiplication denoted as o. Elements of the 
field will be denoted by Greek letters a, /3, 7 . . . £ C, elements of A4 will be denoted by Latin 
letters. Also we assume that A4 contains a unity element which we denote Id. 

All basic objects such as symmetries, first integrals, recursion and Hamiltonian operators can 
be naturally defined for equations on associative algebras. We formulate and illustrate definitions 
for equations on algebra A4 generated by two elements M (variable) and C (constant), but all 
these definitions can be extended to any associative algebra with a finite number of generators 
straightforwardly. 



Definition 1 For a differential equation Mt = F on algebra M we say that G a M is a generator 
of an infinitesimal symmetry if the flows Mr — G and Mt — F commute (i.e if ^ evaluating 
according the equation and ^ evaluating according Mr = G give the same element of M ). 

Remark. In terms of the (non-commutative) differential algebra, equation Mt = F is notliing but 
a derivation Dp of algebra Ai, which maps M to F and symmetry Mr = G is another derivation 
Dg, commuting with Dp- 

To introduce the concept of first integrals, we need an analog of trace, which is not yet defined 
in algebra M. As the matter of fact, in our calculations we use only two property of the trace, 
namely linearity and a possibility to perform cyclic permutations in monomials. Therefore, instead 
of a trace, we define an equivalence relation for elements of in a standard way. 

Definition 2 Two elements fi and f2 of M will be called equivalent, denoted fi ~ /2, iff fi can 
be obtained from /2 by cyclic permutations in its monomials. 

For example, aM oC oM ~ aC oM o M, and of course, any commutator is equivalent to zero. 

Definition 3 Element h of M will be called a first integral of our differential equation, if ^ ^ 0. 
First integrals hi and h2 are said to be equivalent if hi — h2 ^ 0. 

Here there is an obvious similarity with a definition for conserved densities in the theory 
of evolutionary PDEs [^]. In both cases, first integrals and conserved densities are defined as 
elements of equivalence classes, the difference is in the choice of the equivalence relation. 

According to our definition, Rn = M" , n = 1, 2, . . . are nontrivial first integrals of equation 
(H). Indeed, 7?„ 7^ and ^M" = Y2=i^''~^ °C - C o M'^)o M'^-'^ . It is easy to see that 

each element of the sum is equivalent to zero. 

The r.h.s. of equation (H) is a (double) homogeneous polynomial with respect to scaling 
(T : G — > /iG , M — » vM. For such type homogeneous polynomials the exponents of ^ and 1/ we 
shall call weights. As a linear space, algebra M can be decomposed in a direct sum M — (BA4nm 
where a € Mnm if a : a ^ /i^A^a. It is easy to see that 

Lemma 1 Let I be a first integral and G a generator of a symmetry for a (double) homogeneous 
equation. Then any (double) homogeneous component of I or G (i.e. a projection on Mnm) is 
also a first integral or a generator of a symmetry respectively. 

A direct computation of first integrals for ^ with small weights is a simple problem, since 
we have only a few coefficients to determine. For example a general homogeneous polynomial of 
weights (2, 2) is equivalent to 

P2,2 = aG'^ o M^ + PC o M oC o M . 

It is easy to check that — ~ iff a = 2/9 and therefore 

l2,2 = 2G^ o + G o M o G o M (9) 

is a first integral of equation (^). 

In a similar way one can easily find that apart from the obvious orbit integrals Io,n ~ M", 
the following polynomials In.i = G" o M , /i,„ = G o M" and 

I3,2^G^ oM^ + G^ oMoGoM /2,3 = G^oM^ + GoMoGoA^^ (10) 

are nontrivial first integrals of equation (^). It is interesting to note that the space of first integrals 
in each homogeneous component Mkn is one dimensional. We have found the above integrals (^ 
and ( p^ by a simple straightforward computation. We did not use the known Lax representation 
for equation (^) ( not known for all equations), but if we do so, we would find that /„,m is 
equivalent to the projection of {M + C) ("+'") on Mnm- 

A similar computation gives us generators of homogeneous symmetries Sn,™ £ Mnm for 
For instance, to obtain 52.2 we start with a general ansatz 

Q2,2 = aiG^oM^ + a2GoMoGoM + a3CoM^oG + a4MoC^oM + a5MoGoMoG + a6M^oG^. 



It is easy to verify that Q2,2 gives a symmetry for (m iff Q3 = 0:4 = and ai — a2 — —as — — ag. 



Having first integrals and symmetries of equation (|6|) we can use them for the integration of 
Euler's equation. Let C — J^, with J being a diagonal matrix with positive eigenvalues, and 
let M be a skew-symmetric real matrix. Then one can show that H„^rn = trace(7„,m) is a first 
integral of the Euler equation. Of course, only half of the integrals survive such a reduction, since 
Af is a skew symmetric matrix and therefore i?n, 2fc-i ~ 0. Also, the first integrals may become 
functionally dependent. For instance, in the case of the three dimensional Euler top only two first 
integrals are independent. We can choose Ho,2 and //i,2 as a basic set of functionally independent 
first integrals. In the four dimensional case as a basic set we can choose //o,2, H\.2,H2,2 and iifo,4. 
In both cases the Hamiltonian ^ of the Euler equation does not belong to the basic set (it is a 
function of the basic integrals). 

Summarising the above exercise we would like to say the following: The Euler top on an 
abstract associative algebra is not an evolutionary equation and therefore it is difficult, if possible 
at all, to study its symmetries and first integrals directly. Meanwhile, it has an evolutionary 
symmetry (^) which is naturally defined on a free algebra. It is very easy to find first integrals 
and symmetries for (^) (only linear algebraic equations to be solved) . In a standard matrix case, 
symmetries and first integrals of ^ are also symmetries and first integrals of the Euler equation, 
they are sufficient for integration of the Euler equation in quadratures. 



1.2 Other examples of integrable evolutionary equations on M.. 

Equation ^ possesses sequences of first integrals and symmetries. These properties enable us 
to integrate its matrix (finite dimensional) reductions. A natural question is: whether there 
exist other equations which possess such properties? In this section we show that the answer 
is affirmative. The definitions given above proved to be efficient and suitable for handling such 
problems. 

Let us inquire when a general double homogeneous evolutionary equation on M: 

Mt=aM^oC + l3MoCoM + -fCoM^ (11) 

possesses symmetries? According to Lemma 1, we can study symmetries in each homogeneous 
component independently. General equation (^l|) has an obvious sequence of symmetries Si,n = 
Af o C" — C" o M ,n = 1,2,.... In the matrix case these symmetries would correspond to 
the invariance of our equation with respect to similarity transformations commuting with C. 
Symmetries 5*1, „ we shall call trivial and omit them from the further consideration. 

We checked the existence of homogeneous symmetries Sn,m £ Mn,m with n + m < 10 for 
equation (|ll]) by straightforward computations. To do so, we wrote a general homogeneous 
polynomial P„,m £ Mn,m with indetermined coefficients, and commuted equation (|l^) with flow 
Afr = Pn.m- Then we equated the result of the commutation to zero and solved the resulting 
system of algebraic equations for the coefficients of Pn,m and {q,/3,7}. This system is nonlinear 
(all equations are quadratic), but quite simple. All steps in this study are algorithmic. In order 
to perform actual computations we wrote a symbolic code for Mathematica. 



Result 1 Up to scaling of t, every equation of the form possessing a nontrivial symmetry 
Sn,m G M ,n + m < IQ coincides with one of the following equations: 

Mt = [M^,C], (12) 

Mt = [M,MoC], (13) 

Mt = [CoM,M], (14) 

Aft = MoCoM. (15) 

The first equation of the list coincides with (^). Equations ( [l3| ) and ( [l4[ ) are equivalent in the 
following sense. Let us define a formal involution ★ on by 

M* = M, C*=C, (AoB)* = B*oA\ {aA+ 133)* = aA* + (33* , A,3eM, a,(3e<C. (16) 

If an equation Mt = P has a symmetry M-r = S then equation Mt = P* has symmetry Afr = S* . 
The r.h.s. of equation ( p^ is a result of the involution ★ applied to (p^. Similarly, all results 
related to (|l^) can be easily rearranged for ( P^ - 

The form of symmetries for equations (|l3]), (p^ that we found with the help of our computer 
programme has given us a hint for the following statement, which can be easily proven. 



Proposition 1 Equation has commuting symmetries 

A^r„,„ = [M, M"-' o C™] , n,m gN. (17) 

Equation (|7^ /las commuting symmetries 

M^^^ = M oC"" o M , m€N. (18) 

Another interesting problem is to describe equations of the form ( |ll| ) which possess first 
integrals. Straightforward computations lead to the following 



Result 2 Equation ]11\ ) has first integrals of the form In,o ~ M" iff a + /3 + ^ — 0. If equation 
has other first integrals In.m, n + m < 10, then up to scaling of t it coincides with 

Equation (^|) does not have first integrals, and equation ( |l^ ) has only obvious orbit ir iteg rals. 
but both equations have a hierarchy of symmetries. It indicates the "C-integrability" of (y_5|) and 
( |l^ ) (a similar situation for PDEs is described in ^). 

Indeed, if M, C are n x n matrices, then a solution of the Cauchy problem (M(0) — Mo) for 
equation ( p^ is [^: 

M{t) = {e ~ MoCt)-'^ Mo , 

where e is the unit matrix. 

Equation ( [l^ ) in the matrix case can be integrated in the following way. Let Mo be any 
constant matrix and T be a solution of the following linear differential equation with constant 
coefficients 

Tt = MoTC, T\t^o = e (19) 

then M{t) = T~^MoT satisfies equation and M(0) = Mo- Moreover, if T{t,Tn,m) is a 
simultaneous fundamental solution for equation (|l^) and Tr„ „ = MqTC"^ (these equations are 
obviously compatible), then M{t,Tn,m) = T'^MoT satisfies equation A'/t„,„ = [M, M""^C""] as 
well. 

Equations ( |l^ ) and (jl^) have cubic (in M) symmetries. The next natural question is: whether 
there exist other cubic double homogeneous equations 

Mt = aM^ o C + /?M^ o C o M + 7Af o C o <5C o (20) 

which possess symmetries or first integrals. The answer can be formulated as follows. 



Result 3 Up to scaling of t and formal involution every equation of the form ^2C\ ) possessing 
a nontrivial symmetry Sn,m G , n + m < 10 coincides with one of the following equations: 

Mt = [M\ C] , (21) 

Mt = [M, CoM^], (22) 

Mt = [M oC o M, M], (23) 

Mt = M^oC-M'^oCoM + MoCoM'^ ~CoM^,. (24) 

The first two equations (^l|) and ( ^ ) we expected, since they are symmetries of equations ( [l2| ) 
and dl^ ) respectively. 

Equation ( |2^ ) is new, and has symmetries with weights {2k — 1, 1), but does not possess 
symmetries of weights {2k, 1). It is interesting to note that element — satisfies the equation 
Nt = [A^,C] (p^. In the scalar case we would consider these equations to be point-equivalent. 
In the matrix case if we know solutions of (^) we know solutions of equation (^2|) as well, but 
not vice- versa. 

Moreover, we can find an equation for the "n-th root" of N. Indeed, it is easy to show that 
if M satisfies equation 

Mt = M"+^ oC - M" oC o M + M oC o M"^ -Co M"+\ . (25) 

then N — M" satisfies equation (p^. Equation ( |2^ ) has symmetries with weights (fcn + 1, 1) , fc = 
0, 1, . . ., but does not have symmetries of weights {p, 1) where p ^ l(mod(n)). 

In the above list, equation (^^ is maybe the most interesting one. It has symmetries of 
weights (n, 1) for any n. This equation also can be related on a formal level with an extension 
of the hierarchy of equation (p^. If algebra M contains (or is extended by) the inverse element 



N — M ^ , then it follows from (|2^) that Nt — [N ^ , C] and it is a symmetry of equation 
Nr = [^'^C] (Q. Generally speaking, the inverse element may not exist (for example it does 
not exist in the 3x3 matrix case when M is skew-symmetric). 
In a similar way one can verify that equation 

Mt = [M, M" o C o M"] (26) 

is related with (^^ by the substitution A*' = M". Similar to (psj), this equation has lacunas of 
length n in the sequence of symmetries, namely it has symmetries of weights {kn + 1, 1) , fc — 
2,3,.... 



2 Hamiltonian and recursion operators on A4. 
2.1 Basic definitions 

In this section we give definitions for some basic objects such as the Frechet and variational 
derivatives, Hamiltonian operators, etc. suitable for equations on free algebra M. We shall use 
equation (H) to illustrate our definitions. 

For any a G we define operators of left (La) and right (Ra) multiplications which map M 
into itself according the following rule: let b be any element of M, then 

La{b) — aob, Ra{b) — b o a . 

It follows from the above definition and associativity of algebra Ai that RaLt = L^Ra and 

Laob = LaLt, Raob = RbRa, Lca+/3b ~ CtLa + f3Lb, Raa + pb = oRa + fiRb- (27) 

Let us define an algebra O with generators Lm, Lc,Rm, Rc satisfying the following relations: 
RmLm ~ LmRm, RmLc—LcRm, RcLm — LmRc, RcLc = LcRc ■ 

We shall call O as the algebra of local operators. Due to (|2^) any operator of multiplication on 
elements of M can be represented by a corresponding element of O. Let us denote by id the unity 
element of O, i.e. the operator of multiplication by Id. 

The gradation of A4 induces a gradation structure O = ©Ci.m on C If we scale C /iC and 
M vM, then elements of On,m gain the multiplier fi^i/"^. And, of course, if A £ On,m, then 
A : Aip,q — + Aip+n,q+m- wc shall Call elements of On,m homogeneous operators of the weight 
(n, m). 

The Frechet derivative a* for any element a £ A4 belongs to O. 

Definition 4 Let a = a(Af, C) be any element of M. Then the Frechet derivative a* £ O of a is 
uniquely defined by: 

^a{M + eSM, C)|,=o = a,{SM) . 
de 

For example, following this definition we can calculate the Frechet derivative K. of element 
F{M, C)=M^ oC~Co which is the r.h.s. of equation (|): 

-^{F{M + e5M), C))|,=o ^ Y^iiM + eSMf oC -C o (M + eSMf)U^o = 
&MoMoC + Mo5MoC-Co5MoM-CoMo5M = F,{SM) 

Therefore 

K = RcRm + LmRc - LcRm - LcLm (28) 

Let Dp and Dq be two derivations of algebra M corresponding to elements F,G £ M re- 
spectively (cf. Remark after the Definition 1). Then their commutator is also a derivation of the 
algebra A4 

Dk = [Dp, Da] , 

with 

K = G,{F)~F,{G). 



In particular, G is a generator of a symmetry for equation Mt = -F iff 

G,{F)=F4G). (29) 

Usually equation ( p£| ) is called the defining equation for symmetries (c.f. Q]). 
Any derivation Dp of M induces a derivation of O: 

Dp La = La^(F), DpRa = Rat{F) (30) 

for any a £ M. 

The formal involution * induces an involution of O: 

LI = Ra (31) 

It follows from ( ^ ) that i?* = La and in particular _RJ,/ = Lm, Lq — Rc- Operator Q G O is 
called symmetric or skew-symmetric if Q* = Q or Q* — —Q respectively. For example operator 
adM = Lm — Rm is skew-symmetric. It is useful to remember that for any a,h £ M and Q £ M 

aoQ{b) ~ Q*(a) o b. 

2.2 Hamiltonian structures on A^. 

The general Hamiltonian equation on Ai has the form 

Mt = e(grad„(ff(Af,C))) (32) 

where H{M, C) is a Hamiltonian of the equation, O is a Hamiltonian operator. Here we shall 
study local Hamiltonian operators, i.e. assume that O £ C 

For any a G we define the gradient gradj^j(a) £ as follows: 

Definition 5 Let a{M,C),5M £ M. Then grad^j(a[M,C)) is uniquely defined by: 
^a{M + eSM, C%^o ~ SM o gradj^.f{a{M, C)) . 

It is easy to check that if Hi ~ H2 ,Hi £ M then gTa,dj^j{Hi) = gradj,^(_ff2) and in particular 
if Hi ~ then gradjyj(_ffi) = 0. Moreover, the following analom of the well known Theorem (cf. 
§) holds. 

Proposition 2 Let a G A4, then a ~ const, const £ M iff gradj^j{a) = 0. 

For example, let us verify that equation ^ can be written in a Hamiltonian form (^) with 
Hamiltonian H{M, C) = C o and Hamiltonian operator Q = adM (©(a) = M o a -~ ao AI for 
any a G M). It is enough to show that grad^^ (i/(M, C)) = M oC + C o M. Indeed, 

H{M + e SM, C) = H{M, C) + eC o SM o M + eC o M o SM + O(e^) 

and 

^H{M+e5M,C%^o = CoSMoM+CoMoSM ~ SMo{MoC+CoM) = SM ogTa.dj,f{H{M,C)) . 

Definition 6 We shall call O £ O a Hamiltonian operator, if the Poisson bracket 

{a, 6} = gradj^ja o Q^grad^jb), a,b £ A4 

satisfies conditions 

{a,fo} + {6,a} ~ (33) 
{a, {b, c}} + {6, {c, a}} + {c, {a, b}} ~ (34) 

for any elements a,b,c £ M. 



^The above defined gradM is an analog of the variational derivative in the theory of PDEs. 



With the help of the substitution principle g one can show that (^3|) imphes 

e* = -e (35) 

(i.e. B is a skew-symmetric operator with respect to the involution (|3l|)). Moreover, the Jacoby 
identity (^) is equivalent to the following condition (same as (7.11) on page 428 in |^): 

Q2 o XaiQi) + Q3 o XiiQi) + Qi o XiiQa) ~ 0, (36) 

where operators Xi £ O are defined as follows 

X. = De(Q.)(e)- 

It is easy to check that B = ad^t — L/jk — R^k is a Hamiltonian operator for any k. Indeed, 
it is a skew-symmetric operator and satisfies ( ^6| ) since Xi = 0. Operator Oi = adM = Lm — Rm 
also satisfies conditions (js^), (|36|), and therefore, is a Hamiltonian operator on M. Hamiltonian 
operators O and Bi are homogeneous and have weights (fc,0) and (0,1) respectively. For any 
A G M, operator 

B(A) = Bi + A6 (37) 

is a Hamiltonian operator as well. In other words, operators B and Bi are compatible and form 
a Hamiltonian pencil B(A). As usual, the Hamiltonian pencil can be obtained by a simple shift: 
if we replace M by M + XC*" in Bi, we find 

We can employ the method of indetermined coefficients for searching Hamiltonian operators of 
relatively low weights. Conditions (^) and ( p^ ) yield systems of linear and quadratic equations 
for the coefficients of B. 

Result 4 Up to scaling of M and C, every homogeneous Hamiltonian operator of weights (0, m), 
m < 8 and (1, m), m < 7 coincides with one of the following: 

Bo = Rc-Lc, (38) 

Bi = Rm - Lm , (39) 

Ba = LcRm - LmRc , (40) 

B3 = lIiRm-LmRIi, (41) 

B4 = LmLc'LmRm — LmRmRcRm ■ (42) 

The shifts C ^ C + Xld and M —^M + jj,Id in B2 gives a Hamiltonian pencil B2 + ABi + /iBo 
for the Euler top and its hierarchy (i.e. equation (^), for example). Thus equation (^ is a three 
Hamiltonian system. 

A similar shift C — > C + A /d of B4 gives a pencil B4 + AB3. If we assume the element M to 
be invertible, then the following formal change of variables M — relates operators Bi and 

B2 with B3 and B4 respectively. 



2.3 Recursion operators. 

Recursion operators, if they exist, give a convenient way for generating a hierarchy of symmetries 
(see, for instance, [^). For PDEs, most of known recursion operators are pseudo-differential. 
Existence of local, i.e. differential, recursion operators usually indicates on C-integrability of the 
equation j^. Here we begin with a definition of local recursion operators for ODEs on M and 
later extend it to a "non local" case. 

Definition 7 We say that A £ O is a recursion operator for an evolutionary equation 

Mt = F, Me M, (43) 

if it satisfies the following equation: 

At = F.A - AF, , (44) 



where is the Frechet derivative of F 



Operators satisfying equation (Q) act on the space of symmetries of equation (|43[). 

Proposition 3 If A is a recursion operator and G is the generator of a symmetry Mt = G for 
equation then G\ — A(G) ts also a generator of a symmetry Mt^ ~ Gi for equation ^4^- 

Proof. Indeed, let us dieck tfiat (A/t)^ and (Afn)* represent the same element of A4. We have: 
(Mt),, = F^A{G) and 

{Mr,)t = At(G) +A(Gt) = {F,A- AF.)(G) + AG,F = A(G*F - F.G) + F.A . 

Now, due to the condition (p9[), we have (A/t),, — {Mr-^ jt- O 

If we have two operators Ai and A2 satisfying equation then any linear combination 

qAi + /3A2 with constant coefficients q,/3 £ C, a composition A1A2 and, in particular, any power 
A? also satisfy @. 

Applying the sequence A*', fc = 1, 2, . . . to a single symmetry (or to equation ( |43| ) itself) one 
could generate a hierarchy of symmetries. In the case of integrable one-component evolution 
PDEs the whole hierarchy of symmetries is generated by a single recursion operator. In our case 
of equations on Ai the situation is rather different. As we have seen above (cf. (|l^)), integrable 
equations may have two-index hierarchy of symmetries. It gives us a hint that we should have 
two independent recursion operators, and each of the operators raises the corresponding index. 
For example, A = Lm is a local recursion operator for equation ([l3|). Indeed, according (|l^. 



A — Fjvf^c ~ Lmcm- The Frechet derivative of the r.h.s. in this case is 



F, = Rmc + LmRc — RcM — Lmc 

and therefore 

F,A— AF* = {Rmc+LmRc—Rcm—Lmc)Lm—Lm{Rmc+LmRc—Rcm—L 

1+2 



The following sequence IVU, = Lm'^ {M^C - MOM) = M"+^C - M"+'^GM, n e N is an infi 



nite hierarchy of symmetries for equation (p^. It is not the whole set of symmetries (Proposition 
^, (p^)) and we shall see that a "non local" recursion operator is responsible for the remaining 
part of the hierarchy. 

In the theory of integrable PDEs a non local recursion operator, by definition, is a ratio of two 
local (i.e differential) operators. In our case we define the corresponding object as a ratio of two 
operators from O. 

One of possible ways to introduce non local recursion operators for multi-Hamiltonian equa- 
tions belongs to Magri |]lo|]: 

Proposition 4 If Q and Oi are two Hamiltonian operators for the same equation, then their 
ratio 

A = 016"' (45) 

IS a recursion operator. 

Proof. Indeed, any Hamiltonian operator satisfies equation (cf. (7.38), page 449, [^) 

et = F.e + eF.*. (46) 

If Q and Oi are two solutions of equation (^), then it is easy to check that the ratio A — OiO^^ 
satisfies equation (^^. □ 

In order to illustrate this Proposition let us consider equationjq). It has three Hamiltonian 
operators Oo = Lc — Rc = adc (p8|), Oi = Lm — Rm = adM (^) and 02 = LmRc — LcRm 
die]). Consequently, we have two independent recursion operators for (^: 

Ai = {LmRc - LcRM)adll , A2 = {LmRc - LcRM)ad^^ . (47) 

Here we propose a simple generalisation of Proposition and make it suitable even for non- 
Hamiltonian equations. 

Proposition 5 Let Qi,Q2 £ O be two solutions of operator equation 

Qt^F*Q + QP, (48) 

where P £ O is a given operator. Then 

A = Q2Qr' (49) 

is a recursion operator for equation (^j. 



Proof. Indeed, 



At = ^Qr'-Q20r'^Qr' = {F,Q2 + Q2P)QT^ -Q2QT^{F,Qi + QiP)QT'^ = F,K-KF, . 

□ 

E qua tion ( ^8| ) is a generalisation of two well known operator equations. Namely, if P = — F», 
then (^) coincides with equation ( ^ ) for a recursion operator. If P = P,*, then skew-symmetric 
solutions of (^sj) give us Hamiltonian operators (^). Here we would like to emphasise that in 
Lemma ^ we do not specify the nature of the operator P, the only requirement is an existence 
of two (or more) solutions Qi and Q2- Such a construction may occur to be useful for theory of 
integrable PDEs as well. 

It is clear that for a homogeneous equation with P £ Mn.m, operator P must belong 
to On,m-\ and we could try to find simultaneously operators P,Q\ and Q2 by the method of 
indetermined coefficients. Actually, following this simple idea we have found recursion operators 
for equations (p^ ,(|l5|) and (p^). 

Equation ([ISp is not a Hamiltonian system, but it has a hierarchy of symmetries. It is easy to 
find that 

Qn = LmLc , n = 0, 1, . . . 
are solutions of equation (^) with P = RmRc- Therefore 

A = LmLc{Lm)-^ (50) 

is a non local recursion operator for (p^. 

In a similar way one can show that besides of a local recursion operator Ai = La/, equation 
( p^ ) has a non local recursion operator A2 = adMRcad'^l which corresponds to Qi = adM, Q2 ~ 
aduiRc and P = Lmc — LmRc 

It is easy to check that (^), 04 ( ^ ) are Hamiltonian operators for equation Mt = [AI o 
CoM,M] (H). Therefore 

Ai = O463 ^ = {LmLcLmRm — LmRmRcRm){L\iRm ~ LmRm) ^ = (LmLc — RMRc)adj^l 

_ _ (51) 

is a recursion operator for (bq). In the last part of the equality (|5l|), i.e. after a cancellation, we 
have a ratio of Q\ = LmLc — RmRc and Q2 = adM- Neither of these two factors is a Hamiltonian 
operator for (p^), but they both satisfy the same equation (^^, where P = LmLcRm — Lm RmRc 
is a local operator. Moreover, Qs = LmRmRc — LmLcRm is a third independent solution of 
(^). Therefore A2 = Q^ad'^^ is another recursion operator for equation ([23[). 

In spite of a non local recursion operator is formally defined by expression (^9|) in Proposition 
^, its action on symmetries is not well defined yet. Here there are two problems; 

• Action of non local operator (|4s| ) can be correctly defined only on elements of the image 
space AIqi C of operator Qi. 

• If Qi has a non-trivial kernel, then action of A is not uniquely defined. 

Similar picture we had in the theory of integrable PDEs |^], where often there was a problem to 
define the action of D^^ on differential polynomials. 

As a first example we consider an action of recursion operator ( |5o| ) on symmetries of equation 
dl^). In this case the kernel space of operator Q\ = Lm is trivial. The action of (and 
therefore of A) is correctly defined on subset Mq^ = {Mo a; a G M} by Q^^{Moa) = a. Powers 
A*" = LmC'LTj of the recursion operator A are also correctly defined on Mq^ and the r.h.s. 
of equation (113) belongs to A^Qj . Therefore we can generate the following infinite hierarchy of 
symmetries 

S'„,2 = A""^(M oCo M) =MoC" oM, n = 0, 1, 2, . . . 

for equation ([l5|). Equation ( p^ ) has other symmetries, namely the trivial symmetry M^g — 
M o C - C o M and the scaling symmetry M^^ = M + tM o C o M. Action of A (|Ho|) on the 
generator of the trivial symmetry is not correctly defined (i.e. the result does not belong to A^). 
On the contrary, the action of A*" on the generator of the scaling symmetry is correctly defined 
and give the following hierarchy of time dependent symmetries 

Ms^ = A'"^{M + tM oCo M) = Mo C''~^ + tM oC'' o M , fc = 0, 1, 2, . . . . 



As usual (c.f. [^), time independent part of the above hierarchy defines a hierarchy of master 
symmetries. Namely, 

Sm + n,2 ~ Sm,2*{Kn) — K„t{Sm,2) , 

where K„ — Mo C" are generators of master symmetries. 

In the case when operator Qi (see Proposition |^ has a nontrivial kernel we have to impose 
some extra conditions in order to make the action of uniquely defined. As an example, let us 
consider equation with recursion operator Ai (^^. We can choose M*^ , fc = 1, 2, . . . as a basis 
in the kernel space of adM- Applying the recursion operator Ai to the trivial symmetry generated 
by Si,i = MoC-CoMwe obtain 

G = Ai(5i,i) =e2(C)+^afce2(M'=) 

= M oC^ -C^ o + ^ afc(M'=+^ oC -Co M^+^) . 

We see that accounting elements of the kernel space we add symmetries Si^k = oC — C o £ 
Mi^k to G. The function 5*1,1, operators Oi — adM, 62 = LmRc — LqRm are homogeneous and 
if we request that the function G — 5*2,1 ~ Ai(5i,i) to be homogeneous too, we have to choose 
coefficients Uk — 0. It is easy to see that 52, 1 is an element from the image space of adM, and 
therefore 53,1 = Ai(52,i) ~ Af(5i,i) is well defined and is the next member in the hierarchy of 
symmetries of equation etc. Moreover, it is easy to check, that elements Sk,i belong to the 
image of another Hamiltonian operator Qo = adc and therefore the action of operator A2 ( |47| ) 
on these generators is also correctly defined. The whole two-index hierarchy of symmetries for 
equation (|q) can be obtained by the action of recursion operators An,m = A2 A™ on the trivial 
symmetry 

5n,m = Am-l,n-l (5l,l) . (52) 

The same hierarchy of symmetries can be obtained using the first integrals In,m 

S„,m = adMgradj^j{In,m) = 62gradjvj(/„_i,m) . 

In a similar way one can check that pairs of recursion operators for equations ( [l^ ) and ( [23[ ) 
give corresponding two-index hierarchies of symmetries. 



3 Multicomponent equations 

3.1 Breeding and reducing equations 

In the previous sections we were dealing with evolutionary equations on a free algebra M with 
one constant (C) and one variable (Af) generators. Our result indicates (see Result 1) that there 
are only three basic hierarchies of equations and one could think that it may not be worth to 
develop a sophisticated theory to serve these three exceptional cases. In this section we are going 
to demonstrate that many important integrable equations are nothing but particular cases of these 
three equations. 

A way to breed equations is to regard M and C as N x N or even infinite dimensional matrices 
whose entries are generators of a free algebra A and to reduce the system obtained by imposing 
linear constrains on the entries compatible with the dynamics. 

One example which gives the (non-abelian) Volterra equation (^) has been considered in Intro- 
duction. Deriving equation (|) from (^) we have made two steps: we have replaced M, C hy N xN 
matrices whose entries belong to a free algebra and then we have imposed constrains by setting 
some of the entries equal to zero. These constrains are compatible with the dynamics. Indeed, let 
C be given and be of the form (^) . If the initial conditions are such that matrix M has the form (Q) 
then it will remain of the same form for any t. The constrains imposed are not compatible with all 
the symmetries Sn,m of equation (^. Only symmetries 5n,n+i survive under the reduction. First 
integrals of equation (^) can be obtained from the first integrals In,m of (^ by taking a formal 
trace of the matrix corresponding to /n,m (i.e. summing up diagonal elements). For example 

pi = Tr(/i,i) = Tr(AfoC) = ^ Jn o m„ , 



p2 = Tr(/2,2) = Tr(2Af^ oC^ +M0C0 A/ oC) = 

= Jn O Un O J„ O U„ + 2J„ O U„ O U„+l O Jn+1 . 

It is easy to see that Tr(/,i.m) = if n ^ mmod(A'^). 
The modified Volterra equation 

= Uk o {Jk_i o Uk-i - Uk+i o Jk) o Uk, keZpf. (53) 

is a reduction of equation (§3|), where we have to assume Al to be of the form (Q) and matrix C is 
of the form [C]p,q = 5p^2,qJq (here S^^^ = 1 if i = j mod(A^) and otherwise). Again, symmetries 
and first integrals of equation (^sj) can be easily obtained from the corresponding symmetries and 
first integrals of equation (p3|). 

Cyclic reduction (^ is also compatible with equation (p^). The corresponding C-integrable 
system of equations for variables Uk is: 

^Mfc = Mfc o Uk+1 o Jfc+i — Uk o Jk o Uk, k e Zm ■ (54) 

Let us consider the simplest nontrivial case N — 2 and assume that Ji = J2 = Id. In this case 
equations (j54|) and (^) are reduced to 

Ut = uou~uov, Vt~vov — vou, (55) 

and 

Ut = uov — vou, Vt~vou — uov, (56) 

respectively. One more two component equation can be obtained from (^) if we assume — 
3, Ji = J2 ~ J'i ~ Id and uz = —u\ — U2: 

Ut=uou + uov + vou, Vt ~ —V ov — uov — vou. (57) 

All these evolutionary equations are defined on a free algebra A with generators u, v over the 
field C. Elements u and v satisfy a system of equations of the form 

ut — P[u, v) , Vt = Q{u, v) , (58) 

where P{u,v) and Q{u,v) are elements of A. The involution ★ (^^ on A is defined by 

u* = u , V* = V , (a o b)* = b* o a* , a,b £ A. (59) 

Equations, which are related to each other by linear transformations of the form 

u = au + Pv , V — ^u + 5v , aS — f3'y 7^ (60) 



and involution (|59|), we shall call equivalent. 
Equations which are equivalent to 

ut = P{u,v), vt=Q{v) (61) 

we shall call triangular. 

For example, (pq) is a triangular equation. Indeed, after the following change of variables 
u = u,v — u + V we obtain Ut — [ii,!)], vt ~ ■ It is easy to check that equation (pTI) is equivalent 
to 

Ut = v'^ , Vt = u" (62) 
and not triangular. Equations (^^ and ( |55| ) are not equivalent. 



3.2 The simplest classification problem for equations on algebra 
A 

3.2.1 Quadratic equations with a cubic symmetry 



Equations (|55|), ( |57| ) have infinite hierarchies of symmetries and first integrals. It is interesting to 
answer the question whether do exist other quadratic equations 



Ut = diu o u + a2U o V + a^v a u + a^v o 
Vt = PlU O It + I32U O V -\- P'iV o u + PiV o V 



(63) 



which possess symmetries or first integrals? And, if so, how many classes of inequivalent and 
non-triangular equations do exist? A partial answer to these questions is given by 



Theorem 1 Any non-triangular equation \6!^ possessing a symmetry of the form 

Ur ~ jiu o u o u + 721* o u o V + -ysu o V o u + 74U o u o u + 
J5U o V o V + o u o V + 77?; o V o u + jsv o v o v, 

Vt = 5lU O U O U + S2U O U O V + S3U O V O U + SiV o u o u + 

S5U o V o V + Ssv o u o V + 5tV o V o u + 5sv o V o V, 
IS equivalent to one of the following: 



Ut 




u u — u V, 


Vt 




vov~uov + vou 


Ut 




u V, 


Vt 




V u, 


Ut 




u u — u V, 


Vt 




V V — U V, 


Ut 




V, 


Vt 




vov + uov — vou 


Ut 




u V — V u, 


Vt 




uou + uov — vou 


Ut 




V V, 


Vt 




u u, 



(64) 
(65) 
(66) 
(67) 
(68) 
(69) 



It is a remarkable fact, that a requirement of existence of just one cubic symmetry selects a 
finite list of equations with no free parameters (or more precisely, all possible parameters can be 
removed by linear transformations (|60|)). 

The next natural question is whether equations (|64[)-(|69|) have other symmetries, do they 
possess first integrals, whether the corresponding matrix equations are integrable? We have 
calculated all time independent polynomial homogeneous symmetries and first integrals of low 
orders for equations of the list. The dimensions of linear spaces of symmetries and first integrals 
are compiled in the following table: 



Number of Symmetries 



Number of First Integrals 



Order 



6 7 



6 7 8 



9 


10 


11 


12 














1 


1 


1 


1 
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11 


20 


27 


52 
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2 


1 


3 
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(i) 
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ion 
ion 
ion 




1111 



None of the equations from the list has a linear symmetry. The only symmetry of order two 
is the equation itself. 



Equation ( |6^ has a symmetry of order 3 



=uouov , 

Vt=VOUOV~UOVOV. 

(and that is the reason why it belongs to the hst) but it does not have any higher symmetries nor 
any first integrals (at least of orders less than nine and thirteen respectively). 
Equation ( |65| ) has a hierarchy of symmetries 

)=A"('"°M-('"°'°!""1" V n = 0,l,2,... (71) 



where the corresponding recursion operator A is of the form 

A = 

It is easy to verify that 



Ru — Rv, 
0, Ru — Rv 



Pn ^ {u — v)" 



is a sequence of first integrals for equation (|65[). Equation ( pq ) can be written in a Hamiltonian 
form 



where the Hamiltonian operator O is of the form 

RuRu LuLu, LuLy "4" LuRv L^Ru ~h RuRv 



LuRv LyLu LvRu RvRu^ L^Ly R^R^ 

It is interesting to note that equation (^) is a stationary point of 1 + 1 dimensional integrable 
system 

Ut — Ux +UOV , 
Vt ~ —Vx -\- V o u . 

Equations ( |66| ) and (^^ have different sequences of orders of first integrals - it is an indepen- 
dent indication that these two equations are not equivalent, i.e. cannot be related by invertible 
transformation ( |60| ) and involution (^9|). It is easy to check that 

P2n = {UO v)" 



is a sequence of first integrals for equation (|66|) . Its symmetries also can be written in terms of a 
recursion operator 

= A"f \ , n = 0,l,2,... (72) 

\ vov+uov—vou I 



A"( ) , , = 0,1,2,... (73) 



where 



A = 



LuRv 


, LuRu 


— LuLu 


LvRv ■ 


) Lu Rv 









Rv , L 




1 ( 


Rv , L 



Ru Lv 
0, 

Equation ( |69[ ) is, maybe, the most interesting in the list. It is a Hamiltonian equation with 

V -1, 

and Hamiltonian H — |(u^ — u"^). We have seen in the previous section that ( p9[) i s equivalent to 
([5^) which is a reduction of equation (0). Therefore the Lax representation for (M) is known and 



its symmetries and first integrals can be easily found. A straightforward attempt to find a local 
recursion operator (i.e. to find solution of equation (p4|)) gives 



A = 



LvLv ^LvRv RvRvi LuRv LvLu ~h LvRu RvRu 

Lu — Rui \ I Lv — Rv, ~L/u + Ru 
0, Lv — Rv j \ L-u — Rv, L/u + Ru 



(74) 



Equation itself and all known to us time independent symmetries belong to the kernel of this "re- 
cursion" operator A ([t^. Equation (^^ is homogeneous and therefore it has a scaling symmetry 

Applying A to this scaling symmetry we obtain a cubic symmetry 

j ^ / 11 + tv^ \ f ~2u ouov + 4uovou — 2vouou \ , , 

Vt J \ V + tv? j \ 2uovov — 4vouov + 2vovouJ 

of equation (|69|). 

A systematic description of symmetries and first integrals for equations ( p7[ ) and (^^ is still 
an open problem. 

3.2.2 Quadratic equations with a quartic symmetry 

Integrable equation (^) does not possess a cubic symmetry and therefore does not belong to the 
list (p^-(|69|), but it has a quartic symmetry 

2 2 

Ut = U O V o u — U O V o u , , , 

2 2 

Vt=vouov— vouov. 

Symmetry ( [t^ ) and higher order symmetries of equation ( |55| ) can be generated by the following 
recursion operator: 

LuLv: 

0, LvLu 



A = 



Proposition 6 Apart from ^6q)-^6q) the list of quadratic equations on A which possess a quartic 
symmetry includes the following meqmvalent equations: 

Ut = uou~vou, , , 

Vt=vov — uov; 

m = -uov (78) 
Vt=vov + uov; 

ut^^vou, 
Vt=vov + uov\ 

Ut = uou~uov~2vou, , , 

Vt=vov — Zuov — vou; 

Uf— uou-'2vou, , , 

Vt = V o V — Zv o u ; 

Ut ^uou — 2uov, , , 

Vt=vov + 4:Vou. 

Here we do not claim that the list presented above is complete. Some extra work is required 
to cast Proposition ^ into the form similar to Theorem 1. 

It is interes ting to look at the sequence of dimensions for symmetries and first integrals of 
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Equations (|77|), (|78|) and ([79|) have identical sequences for dimensions of symmetries and first 
integrals. Nevertheless we have checked that they are not equivalent, i.e. cannot be related by 
linear transformations (^) and involution (|59|). 

Equation ( js^ possesses a symmetry of order four but it seems it does not possess any other 
symmetries or first integrals. 



Instead of a conclusion 

In this paper, we have made an attempt to bring the rich hidden structure of integrable PDEs M] to 
a new domain, namely, to differential equations on free associative algebras (see also ||l^, p^, p^]). 
We have formulated basic definitions and shown their correctness and efficiency. Like integrable 
PDEs, ODEs on associative algebras may have infinite hierarchies of symmetries and first in- 
tegrals and that asserts an algebraic definition of integrability. In the case of finite dimensional 
(matrix) representations of the algebra, the corresponding (matrix) systems of ODEs inherit these 
symmetries and first integrals and can be integrated in quadratures. 

This study raises a lot of questions and open entirely new area for research. For instance, 
we foresee that the list of integrable equations on associative algebras with constrains should be 
a lot bigger. Quantum problems, with some commutation relations, naturally fall in this class. 
Another important and promising problem for further research is a systematic study of reductions 
of equations on associative algebras and, in particular corresponding matrix systems of ODEs. 
We expect that infinite dimensional realisations (for instance by operators in a Hilbert space) and 
their reductions may be of interest as well. We would like, also, to look at the theory of classical 
integrable tops from the point of such a componentless approach. That would require to study 
equations on associative algebras with a few constant elements related by algebraic constrains. 
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